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We discuss the prospects of observing double parton scattering (DPS) 
\ processes with purely leptonic final states at the LHC. We first study 

\ same-sign W^W^ pair production, which is particularly suited for study- 

ing momentum and valence number conservation effects, followed by dis- 
cussions on double Drell-Yan and production of J/ip pairs. The effects 
of initial state and intrinsic transverse momentum smearing on pair-wise 
transverse momentum balance characteristic to DPS are studied quanti- 
tatively. We also present a new technique, based on rapidity differences, 
to extract the DPS component from a double J/ip sample recently studied 
><| : at the LHCb. 
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1 Introduction 



Due to the composite nature of hadrons, it is possible to have multiple parton hard- 
scatterings, i.e. events in which two or more distinct hard parton interactions occur 
simultaneously, in a single hadron-hadron collision. For a given invariant mass, such 
cross sections tend to increase with collision energy due to the rapidly increasing par- 
ton fluxes when successively lower momentum fraction x is being probed. The high 
collision energies at the LHC thus provides a valuable opportunity to observe mul- 
tiple parton hard-scatterings. In particular, many double parton scattering (DPS) 
processes involving leptonic final states could become accessible for the first time. 
These include double electroweak processes, for example production of same-sign W 
pairs (VF^VF^) and double Drell-Yan (DDY) interaction, and also pair production 
of J/ip^s. Compared with the DPS processes already observed, namely final states 
involving 4 jets (at the AFS collaboration at the CERN ISR [Ij), and 7 + 3 jets (at 
the CDF [2] and the DO [3] collaborations at the Fermilab Tevatron), properties of 
the leptonic final states can be measured much more precisely. These processes also 
involve different scales and initial state partons, and hence provide complementary 
information on the non-perturbative structure of the proton to the information de- 
rived from other DPS reactions. It is therefore important to study properties and 
prospects for observing various DPS processes in detail. 

The general expression for the DPS cross section cr^/'^^ is given by 

^fyCi) = ~^ J dxidx2dx\dx'2d% 

i,j,k,l 

xTij{xi, X2, b; ti,t2)Tkiix[, x'^, b; ti, t2)crf^{xi, x[)afi{x2, X2) , (1) 

where Tij{xi,X2,b]ti,t2) is the generalised double parton distribution function for 
partons i, j with momentum fractions Xi, X2 at scales ti = ln{Ql), t2 = ln{Q2). The 
two partons are separated by a transverse distance b. The scales ti and ^2 are equal 
to the characteristic scales of subprocesses A and B respectively. The quantity m is 
a symmetry factor that equals 1 if A = B and 2 otherwise. 

For processes that probe small x values, different partons may be expected to 
scatter independently to a good approximation. In this limit, we have 

Tij{xi,X2,b) = Dij{xi,X2)F{b) , (2) 
Dij{xi,X2) = Di{xi)Dj{x2) , (3) 

where the scales are implicitly set to equal values (ti = ^2)- The first expression fac- 
torises Tij into a longitudinal double parton distribution (dPDFs) Dij and a (fiavour- 
independent) transverse distribution F{b). In the second expression, Dij is further 
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factorised into two single parton distribution functions Di and Dj. Using these as- 
sumptions, crf/'^) can be written as 

(4) 

The quantity a^s is expected to be energy and process dependent, and is one of 
the DPS properties that requires a more precise experimental measurementEl Ba- 
sic sum rule constraints, namely momentum and valence number conservation, may 
be included via the longitudinal dPDFs Dij and could be probed in particular pro- 
cesses. In the following, we discuss the prospects of obtaining this information from 
the W^W^, DDY and double J/ip processes at the LHC. The discussion is based 
on Refs. [H El E]. We refer readers to other contributions to these proceedings for 
phenomenological studies involving jets and recent theoretical and experimental de- 
velopments. 
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2 Same— sign pair production 

It was first sugggested in Ref. [7] that same-sign W^W^ process might be a clean 
channel for the observation of DPS. The irreducible background to the process, namely 
single parton scattering (SPS) production of a VF^iy^ pair must be accompanied by 
two additional partons in order to conserve electromagnetic charge. These extra par- 
tons provide excellent handles to separate the DPS signal from the SPS background. 

It was subsequently pointed out in Refs. [6l |8] that in certain regions of final state 
phase space this process could be particularly sensitive to momentum and valence 
number constraints that must operate at some level on the two parton PDFs. The 
relevant region is the one in which both Ws are produced in the same forward direc- 
tion, since this configuration favours extraction of two large x (valence) quarks from 
one proton. Compared with simple factorised models, the sum rule constraints would 
suppress such configurations. The charged lepton r] asymmetry 

^ ^ (^{Vh X Vh < 0) - (y{rih X Vh > 0) 
" cx{7]i, X rji, < 0) + a{7]i, x r]i, > 0) ' 

where rji^ are the lepton pseudo-rapidities, should hence be positive for |?7zJ, |?7i2l > 
^min some minimum pseudo-rapidity cut r;™™, and should also increase with r;™™. 

The functional dependence of a^^ on r^™™ is plotted in Figure [H assuming that 
proton-proton DPS can be described in terms of dPDFs and taking various different 
forms for the dPDFs. The predictions from the GS09 dPDFs [8J in which momentum 

*For concreteness, in the following numerical studies we shall use acs — 14.5 mb, the value 
obtained by the CDF experiment [2] . 
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Figure 1: Charged lepton pseudo-rapidity asymmetry as a function of ?7™™, the mini- 
mum lepton 7] cut, for positively-charged leptons from DPS in pp colhsions at i/i = 14 
TeV, evaluated using different dPDF models. 

Cuts 
\Vi\ < 2.5 

20 <p^T< 60 GeV 

> 20 GeV 
OS lepton veto 

Table 1: Selected cuts (left) on the same-sign muons, and cross sections [fb] (right) 
after cuts for signal and background same-sign W production, including branching 
ratio into muon. 

and valence number constraints are implemented, are compared with those from sim- 
ple "MSTW„" sets, defined as Dij{xi,X2) = Di{xi)Dj{x2)0{l — xi — X2){1 — Xi — X2)^ 
for n = 0,1,2, in which momentum constraints are implemented very crudely and 
valence number constraints are not included at all. It is known that there are theoret- 
ical problems in describing DPS using dPDFs (see e.g. [9j) and none of the dPDF sets 
used take account of the potential contributions to DPS starting from 2 or 3 nonper- 
turbative partons in a correct way. However, such contributions play a subdominant 
role in determining the shape of a^, (^7™™), and the main force shaping this distribution 
is in fact the inclusion of basic momentum and number sum rule constraints, which 
are contained at least in an approximate way in GS09. 

To see how well the DPS W^W^ signal can be extracted from background, we 
perform a parton level study, relevant for the general purpose detectors ATLAS and 
CMS, to look for same-sign muon pairs. We include diboson {W'^Z/'j* , Z/'j*Z/'j*), 
heavy flavour (66) production as well as the irreducible W^W^ + j's SPS background 
discussed above!!] The cuts and resulting cross sections are displayed in Table [H We 

^We refer readers to the original paper [6] for a technical account of the simulation. 
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Cuts DDY cross sections [fb] at LHCb 



1.9 < 7] < 4.9, pt>1 GeV 




DPS 


SPS 


m^+f,- > 4 GeV 


7 TeV 


0.08 


0.43 


9.2 < m^+^- < 10.5 GeV veto 


14 TeV 


0.16 


0.68 



1^4(1 < 40 GeV veto 



Table 2: Cuts (left) and DPS and SPS DDY cross sections [fb] (right) for pp collisions 
at 7 and 14 TeV. 

see that a small excess can be expected, with the background dominated by diboson 
production. There are however additional handles which can help distinguish the 
signal from background. See Ref. [6j for more details. 

3 DDY and double J/i/j at LHCb 

One of the main characteristic in DPS is the so-called pair-wise balancing, in which 
the final states from the two hard scattering processes have zero transverse momentum 
at parton level. This has been used to help identify the presence of double parton 
scattering in previous experiments. If all four DPS final states are charged leptons, 
in principle pair-wise balancing could also be observed. 

In the following, we focus on events with four muon final states, forming two 
opposite sign (OS) muon pairs. We look in the low invariant mass region, where 
the DPS to SPS ratio is expected to be larger. LHCb has excellent low p^ muon 
acceptance, which can go down to ~ 1 GeV, and muon identification in the low 
mass region, making it well suited for studying 4-muon DPS events. In the low mass 
region, the factorisation into two single 2-to-2 hard scatterings in Eq. H] is likely to 
be a good approximation, and so is used in the numerical analysis that follows. 

We first discuss DDY. We use Herwig++ v2 . 4 . 2 [TU] to generate the DPS signal, 
and MADGRAPH v5 . 1 . 2 . 4 [TT] to generate the SPS 4-muon background, which is then 
interfaced to Herwig++ for parton showering. We also include a Gaussian intrinsic 
Pt smearing, parameterised by the parameter a = 2 GeV, to provide a more realistic 
simulation. The cuts and the resulting cross sections are displayed in Table [21 We see 
that while at 7 TeV, the cross section is too low for an expected integrated luminosity 
of 0{1) fb~^, it might be possible to observe DPS events at the 14 TeV LHC at high 
luminosity. 

In principle, the DPS signal and SPS background can further be distinguished 
by pair-wise balancing. However at low invariant mass, the presence of initial state 
radiation (ISR) and intrinsic pt smearing significantly affect the balancing property. 
Due to the ambiguity in grouping the 4 muons into 2 OS pairs, a pair-wise balacing 
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variable (S) 



is used to group the muons into 2 OS pairs by minimising 5*. The 5* distributions 
including different radiation effects are displayed in Figure |2] (left plot). Clearly, 
these distributions depend sensitively on the radiation effects, with the SPS and DPS 
distributions becoming more similar after these effects are included. 
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Figure 2: Distributions of the pair-wise balancing variable S (left) and the rapidity 
difference \ Ay{fi~^fi~ , yU+/i~)| (right) for the DPS and SPS processes including different 
radiation effects. In the figure, PL stands for "parton level", PS stands for "parton 
shower" , and a indicates the inclusion of intrinsic pt smearing. 

On the other hand, longitudinal kinematic variables are expected to be less sen- 
sitive to the radiation effects, which primarily affects kinematic distributions on the 
transverse plane. One such variable is the rapidity difference \Ay{fi~^fi~,fj,'^n~)\ 
between the two OS muon pairs specified by minimising S, which has the addi- 
tional advantage of being invariant under longitudinal boost. The distributions of 
|A?/(yU+/i~, /i+yU")| are displayed in Figure [2] (right plot). We see that both the DPS 
and SPS distributions are much more stable against ISR and intrinsic pt smearing. 
Also, the fraction of DPS events increases with |Ay(/i"'"/i~, /i+/i~)|, making it an ex- 
cellent variable to distinguish DPS from SPS events. 

The above observations can be applied to double J/ip production. Compared 
with DDY, this process benefits from the large single J/ip cross section, while the 
theoretical description of the production is an active area of current research. In 
double J/iIj production there is no ambiguity in grouping the 4 muons into 2 OS 
pairs, as the correct pairing should have /i"*"/!" invariant mass close to the physical 
J/ip mass. The set of cuts and the resulting cross section is displayed in Table [3l 
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Cuts Double J/ip cross sections [pb] at LHCb 



1.9 < r/ < 4.9 




DPS 


SPS 


Pt > 1 GeV 


7 TeV 


3.16 


1.70 




14 TeV 


7.69 


2.62 



Table 3: Cuts (left) and DPS and SPS double J/ip cross sections [pb] (right) including 
BR{J /ip — )■ fi~^fi~) for pp collisions at 7 and 14 TeV. 
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Figure 3: Distributions of invariant mass mj/^j/^ (left) and variation of the cross 
section as a function of minimum rapidity difference \Ay{fi^ fi~ , ^'^ (right) for 
the DPS and SPS processes. Note that slightly different cuts, motivated by the cuts 
used in the recent LHCb measurement [12], have been used to obtain the above 
distributions, c.f. jHIS]. 

In fact, recent results from LHCb [12] might already indicate the presence of 
double J/ip events from DPS. Assuming cXefr = 14.5 mb, the theoretical cross sections 
for the SPS and DPS double J/ip processes are similar. However the mj/^j/^ invariant 
mass distribution is different, with the DPS distribution peaking at slightly higher 
values. In Figure |3] (left plot), we show that combining contributions from DPS and 
SPS might provide a better fit to data. 

The large double J/ip cross section allows extraction of the DPS events from the 
SPS background by imposing a cut on minimal \Ay{fi~^ fi~ , fi~^ In Figure |3] (right 
plot), we show the variation of cross section as a function of min |A?/(/i+/i~, yU"^/i~)|. 
Increasing min \Ay{fi'^fi~,fi^fi~)\ can result in much higher DPS fractions. With 
more upcoming LHC data, this could be an excellent tool in establishing the presence 
of DPS double J/ip events [3]. 

4 Summary 

Measurements undertaken at the LHC will be crucial for improving the theoreti- 
cal description of multiple parton scattering, in particular double parton scattering. 
We have studied the prospects of observing DPS with leptonic final states. Both 
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same-sign and Drell-Yan pair production can be considered standard candle 

processes at hadron colliders, as the leptonic final states provide clean experimen- 
tal signatures, while the theoretical predictions of the 2-to-2 subprocesses are under 
control. The cross sections for these processes are of (9(0.1) fb, and so will require 
high luminosities to obtain unambiguous signals. On the other hand, double J/ if] 
production has a much larger cross section, and a significant DPS component may 
already be present in a recent LHCb study. With more data in the 7 TeV run, the 
DPS component can be disentangled from the SPS component using the rapidity 
separation between the two reconstructed J/ip^s. 
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